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CALCULATION  OF  COSMIC-RAY  INDUCED 
SOFT  UPSETS  AND  SCALING  IN  VLSI  DEVICES 


INTRODUCnON 


Progression  of  VLSI  (Very  Large  Scale  Integratipn),  circuitry  to 
smaller  feature  sizet  substantially  increases  the  i>robabiliiy  of  s^ft  upsa^ 
induced  by  the  penetration  of  energetic  cosmic-ray  particles  thr^igh 
device.'  These  devices  can  change  their  logic  state  withgut  peneNoent 
damage  to  the  device  when  a  densely-ionizing  particle  deposits  a  giiintity  of 
charge  at  a  node  (an  NDS  capacitor,  for  example)  that  is  ebmper^le  to  the 
quantity  of  charge  representing  the  logic  state. 

An  exploratory  study  is  needed  at  this  time  to  estimate  the  effect  of 
scaling  of  microelectronic  devices  to  smaller  sizes  on  expected  ;SOftHjpset 
rates  in  the  cosmic  ray  environments  encountored:  by  satellite^^  fi 
definitive  calculation  requires  knowledge  of  the  ,way  ipiqi^leetrQiHc 
technology  will  proceed  on  scaling  to  smaller  feature  sizes.  However,  some 
limiting  cases  can  be  studied  which  will  Indicate  whether  soft-upsft  rates 
become  impossibly  large  on  further  scaling  or  wither  tolerable 
rates  are  to  be  anticipated  as  we  evolve  to  the  VliSIC  (very  h‘1^  spejd 
Integrated  circuit)  era.  Such  an  exploratory  calculation  h|»  bito  ctf^teilt 
out  by  Burke. z  The  calculations  reported  herein  are  to  extensibir  ^^^^ 
Burke  calculations  in  that  our  calculations  are  based  on  extended  ceamii^^ 
environments  recently  developed  by  the  NRL  Laboratory  ^ 

Physics, 3  and  our  calculations  toplOy  exact  evaluation  bf. 
chord-length  distributions.  Scaling  of  the  devito  sti^tures  is-caj^totolp 
feature  sizes  of  approximately  0.1  micrometer,  arul  venetion  In 
laws  is  explored.  , 

The  study  by  Burke  parametrizes  a  cosmic-ray  LET  spectfito 
Heinrich. A  The  Heinrich  spectrum  only  includes  cosmic-rav  toaienehis  wito 


Heinrich, A  The  Heinrich  spectrum  only  includes  cosmic-raty 
6  ^  Z  <  26,  whereas  the  NRL  spectrum  includes  a  larger  ram 
28.  In  addition,  the  NRL  analysis  accounts  for  tiiM^II 
producing  data  for  three  cosmic  ray  environments: 
solar  minimum  activity,  and  a  90  percent  worst  cdw  ctod 
percent  worst  case  is  the  envlrofiaent  moit  rttooneble  to  toe 
studies  for  sitelllte-borne  eldctrohics  sitod  Ibis  Jetoi 
by  definition,  will  only  be  exced^  10  pifttot  toe  tlm 
soft-upset  rates  have  been  calcdlated  for  m  di|% 
protpnt  in  the  radiation  belts  at  an  aHltudl  Of  iiO„iietofl 
contribution  to  soft-upset  rates  by  prOIgn^-itoucfi 
subject  of  a  continuing  investigation.  Tm  ptostot  ealcuilt 
the  dependence  of  expected  soft-upset  rates  on  environment. 
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Exact  evaluation  of  the  Integral  chord  length  dlstrfhutloiit  |if 
conparlson  with  ^roxlnate  expressions  which  have  been  used  by  o^hsifs  lit 
shown  for  selected  exanples. 

Evaluation  of  soft  upset  rates  for  the  sane  randon  access  nsnorlet 
(RANs)  that  have  been  studied  by  Burke  has  been  per^naed  using  the  seM 
geonetrical  paraneters  and  critical  energies  In  order  to  facHltete 
conparlson. 

II.  EMVIROHMEMTS 

The  LET  spectra  for  the  cosn1c«ray  environnents  were  derjye^.fron  Vn 
work  of  Adans  et  a1.3  We  considered  the  pure  gallctlc  c<^c  rhy 
envlrorment  at  tines  of  ninimun  and  maxlmun  solar  activity.  The  .t^T, spectra 
for  these  two  environnents  are  shown  In  Figures  1  and  Z  respidfl^y,  For 
the  BO  percent  worst  case  spectnait' tie  have  IndudiNy  a  90,|ier<^t  iorifi 
contribution  of  low  energy  particles  ffon  solar  and  intehiilam^y  adtiyfly 
In  addition  to  galactic  cosnic  rays  aS  they  appear  during  ^  per1di  |^ 
mlnlsaai  solar  activity.  Only  10  percent  of  the  tine  shouTt  condltlhj^  In 
the  Interplanetary  medlun  be  nore  hostile  tor  nlcroelectrionli  coipoilnts* 
This  90  percent  worst  case  spectrun,  of  extrene  Inportince  fof  estinatidn  of 
satellite  vulnerability.  Is  shown  In  Figure  3. 

A  copy  of  the  LET  spectrun  for  cosnic  rays  due  to  Heinrich^  is  ^MMh 
for  reference  In  Figure  4.  This  LET  spectruh  Is  very  close  to  the,  splir 
naxinun  siMctrun  generated  at  NRL.  ' 

Figure  5  shows  the  LET  spectrun  for  protons  for,  a  600 
oH>1t  In  the  proton  radiation  belts.  The  proton  SiHi^tjf^  nai 
using  the  radiation  belt  proton  spectrdn  of  APWW  1nf(lir4tw^ 

63®  nn  kn  orbit,  and  allowing  for  the  Sh1e1di)%  of  4  typit^i 
spacecraft.  The  proton  energy  s^triin  was.  cfwvtim  lb  m  ^spbiEm 
using  the  tables  of  Wlllllnson  and  Bbujot.7  thlirspe^ruM  lncilfS^wE 
the  direct  Ionization  effects  of  the  p^tchis;;  MWl  HI 

lolffaatlon  df  any  reaction  produicts.  As  pointed -‘'out:  lb-  reyireb<^.;lLlni^^ 
the  proton  reaction  products  can  easily  m 

devices.  So  that  the  actual  upset  irateb  are  nueh^  hlg^  at  laIttL  iliiK 
factors  than  given  by  predictions  using  this  spectrue.  This 
for  showing  the  worst  case  upset  rates  as  scaled  devices  beCfSse 
the  direct  Ionization.  ' 


'■ ,  -  - 


OlSCUSSIOh  OF  CALCULATION 


Technlgees  orlelnaUy^^  1h 

utilized:  In  cafeuiat^  pf ^ej^eclid  soft  ii 
•redfo^  his  shoun  NOe  td  i4li^  ItUll 
randot  traversals  of  coheex  bodies  by  sm^ 
eapected  tofuei^  ratai.  The  pf«Mlt 
here  closely  felloMi  thit  Of  OerteX  in 
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and  diffused  Junction  test  _  _  _ _ _ _ _ _ 

significantly  less  charge  fron  the  same  ^article  than  dl^it*^ 
devices.  The  energy  deposited,  E(j^,  la  given  by 

i4 

Edep  •  Lpt*  (1) 

where  L  Is  the  linear-energy  transfer  for  the  particle  (LET),  e  1$ 
density  of  the  material,  and  i  Is  the  chord  length  for  the  particular 
traversal  through  the  sensitive  volume.  Thi  nuatber  of  electron-hole  pairs 
created  Is  Ed«p/Ee-h*  where  Ee_h  Is  the  energy  required  to  create  an 
electron-hole  pair,  3.6  eV  In  the  case  of  SI.  hhen  Edeo  sroi^  M 
aE,  a  critical  energy,  a  sufficient  nunber  of  e1ettr(^-4i^  |riinr  W 
to  cause  an  event.  Thus,  an  event  can  occur  If  the  chord  leng^  as^llihd 
with  the  particular  passage  of  a  particle  Is  greater  than  where 

*  min  “  aE/pL.  (2) 

It  Is  shown  by  Kellerer.^Z  that  In  an  Isotropic  uniform  field  of  fluenee 
a,  the  expected  number  of  chords  through  a  convex  body  of  surface  S  Is 
sa/4.  If  this  number  Is  multiplied  by  C(tfBi4)),  the  sum  probability  that 
the  chord-length  through  the  sensitive  volume  Is  greater  t^.. 1 4 
the  number  of  events  that  deposit  an  energy  greater  then  altls  mfliMd. 
For  a  continuous  spectrum,  he  is  obtained  by  Integration  of 
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In  this  Integral,  #(L)  Is  obtalntfd  from  #(E)  using  the  transforietllNi 


The  SMin  dependence  In  C(  convertid  to  a  di|ai 

using  equation  (2)  so  that  tile  1ntefratlM>ovtr  L  cin  bo  piv^ 
Tower  limit  of  Integration  Is  the  rninlsil  spTuo  of  L  thifl  oih 
event,  ^ ^  t k 

L0  •  -.r 

where  Saax  ^  diagonal  of  the  parol 
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U  th®  naxlmum  VAlii*  of  L  Incluiljid  In  #(1).  Jn  ^  (Wfwt . ^ ^ 
1S®»  of  L  (C4l7  it  X*)  cowtafrW  trir  TO  U 
1$  to  b®  *®ro  i®  th®  int#rvtU  to  1,^- 

Ah  jiUerhhtiv®  fohiuUtlQn  to  cafdiitii® 
to  oquailbn  (3)  has  b®®ri  (is®d  h!if  Hclc®T  apd,  Bianwwcl,!*  itiis  tbninwtofi 
uses  ♦4t)t  tn®  thte^al  LET  ^ictruir  th<*  .  ^(  *  >•  the  4WUr®f#M«l 
chord-langth  dlstrlbutloh.  These  (luentftlas  are  related  to  d(L)  aMd  CliKiy 

'  '  »  -■•■'•  '.>i 

/hiiax  , 

#(x)dx.  (7) 

^max 

c(*)  -  f  f{x)dx.  m 

£ 

The  equivalence  of  the  two  formallsais  Is  easily  shown  by  expanding  equation 
(3) 


L  £ 

c  /•  nax  ^x 

Ng-f  y  TO  J  f(f)d£'dL. 


where  s  •  aC 

(  •  -• 

Interchanging  the  order  of  Integration*  with  the  appropriate  change  of 
lleltStwe  have 

’  L 

N®  - 1  /  f(t)  J  #(L  )dL  d£  (10) 

£^  L(£)  x 

Here,  L(£)  «  aE/p£,  the  nlnlanae  value  of  LET  which  will  deposit  tn  the 
sensitive  voIuon  for  that  value  of  L»  and  £q  »  a£/pL|Muc  Using  (7) 
above. 
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Kellerer^Z  does  not  fwmish  explicit  expressions  for  tho 
chord-len^  distributions  thet  appear  in  hli  theoretical  work  on 


HoMver.  Bradford  Jmo 
ColonwM  4«r  the  dirfferwitiel 

tMO«df«en$tone1  rectangle  together  with  tlMitfeciiaHM  ^Ielgiidd 
relatlng^  the  various  two*  and  tht^e-dleanalonal^^choi^^eegth '^litrfttd^^ 
to,  derive  an  aopiroxlmte  expreaslon  for  C(4>,h  the  4«t#^fc^dior*siii^ 
distribution  for  a  rectangular  parallelepipeds  .  lhe,^rox««ptte  1t^c1alild 
to  be  useful  for  b/a  or  c/a  x  1.  This  appi^lawte  eweledilaii  iasf'^jCf  ftjK’is 
very  useful  since  It  furnishes  an  expression  for  4ft)  ie  elosed  fans. 
Burke»  In  the  worked  cited  above*  has  fwsed  a  further  approxllaiyon  t».4^ 
Bradford  expression  which  1$  In  e  oonvenient  form  for  hand  oaleulatfon^lS 
This  approximation  Is 

C(t)  -  0.75  (a/ 1)2*2  i>a  (12a) 

C(i)  -  1  -0.25  (i/a)  i<  a.  (I2b) 

These  approximate  results  are  very  useful,  yielding  values  of  6(t)  gulte 
close  to  the  exact  result  except  that  C(a  )  does  not  go  to  zero  correctly  for 
large  1  near  imax* 

An  exact  expression  for  the  differential  chord-lenfth  distrjbutfgn, 
f(  i),  due  to  Petroff  Is  contained  In  a  papar  by  Pickel  and  BlandfQrd.P 
He  have  evaluated  C(  i)  by  numerical  Integration  of  the  Petroff  result* 
utilizing  equation  (8).  This  exact  result  for  the  geometries  used  In  tiie 
soft-upset  calculations  Is  compared  with  an  evaluation  utilizing  Bradford's 
approximation  (Approximation  1),  and  the  approximation  used  by  Burke 
(Approximation  2).  We  have  also  compared  the  Bradford  approximation  with  an 
exact  calculation  for  a  case  studied  by  Zlegler.^^ 

Comparisons  between  the  exact  calculation  of  C(i)  and  Approxlamtlen^  1 
arc  shown  In  Figures  6  to  9.  He  note  that  results  are  almost  -IdpnIilciJ 
small  1  until  the  first  discontinuity  In  C(i )  Is  reached*  aM  thi^ 
results  are  In  quite  good  agreement  beyond  that  until  i  reaches  <a  yilue 
larger  than  the  last  discontinuity  In  C(i)*  For  these  lar^  valuqS  4f  s  * 
Approximation  1  for  C(i)  does  not  go  to  zero  properly  bt  jUgir*  ,1^ 
examples  shown  In  Figures  6  to  8  satisfy  the  criterion  given  by  liiHipniird'fpr 
the  usefulness  of  the  approximate  formula.  Figure  9  shpws  a  saii^ulirt  the 
parallelepiped  approaches  a  cube*  the  case  studied  by  Utgler.  In  ease 
the  disagreement  at  large  chord  lengths  Is  greater.  r 

Approximation  2  Is  compared  with  the  exact  calculatlea  faf^^s^l.  W 
Figures  10  to  12.  Again*  the  results  using  this  ahif^|i#le9i  4ie 
reasonably  close  to  the  exact  results  except  for  the  longer  chord  tenglhi* 
The  degree  of  agreement  seems  to  ba  depepdent  on  hqm  well  the  oondltlon  b/a 
2l  3  Is  satisfied.  '  , 

IV.  DEVICE  SCALING  AND  CALCULATION  OF  SOFT-UFSEt  RATE  '  ^ 

A.  SCALING 

Soft-upset  rates  (n  the  cosmic-ray  anvlrofsien^  have  been 
calculated  for  the  particular  memory  devices  studied  by  Burke.^  The 
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scaling  scenarios  rep<Mrtcd  by  Burke  are  extimded  t»  smaller  ftature  files 
and  soft-upset  rates  are  calculated  for  ^  vdrlous  envlroiwients  dll^^ 
In  Section  11.  Furtbemere*  where  equation  3  Is  utilized  to  catetlati  the 
number  of  events  that  deposit  energy  greater  then  aE  by  cosmic  ra|e  or  by 
direct  lonizetlon  from  protons  In  the  radiation  belts,  exact  evaluatloh  of 
C(A)  Is  employed.  Mhen  scaling  to  smaller  device  size,  equation  3  can  he 
considered  to  apply  to  the  reference  devices,  4K  memories.  It  Is  not  clear 
at  this  time  how  the  technology  will  evolve  when  VLSI  (tevices  are  scaled 
down  to  smaller  feature  sizes.  One  scaling  scenario  ^at  cm  be  considered 
Is  scaling  according  to  the  model  of  Mead  and  Conway. In  this  model, 
all  of  the  linear  dimensions  are  reduced  by  a  scaling  factor  «, 


a/a 

(13a) 

b/a 

(13b) 

c/a 

(13c) 

Furthermore,  all  voltages  are  scaled  down  by  dividing  by  the  same  scaling 
factor  a,  thus  keeping  all  electric  fields  In  the  device  constmt.  With 

these  assumptions,  the  stored  charge  representing  a  bit  scales  as 

This  follows  from  the  fact  that  the  capacitance,  C,  being  proportional  to  an 
area  divided  by  a  separation  distmce,  scales  as  Since  Q  •  CV,  Q 

scales  as  a*^.  Critical  energy  Is  directly  proportional  to  critical 
charge,  I.e.,  aE  >  tQ  where  c  •  3.6  eV/eh-palr  In  silicon.  Tims  the 

critical  energy  for  upset  scales  as 

aE*  •  aE/a2.  (14) 

The  Mead  and  Conway  scaling  scenario  Is  rather  simplistic  md 
cannot  be  expected  to  apply  Indefinitely  as  devices  get  smaller'  and 
smaller.  For  example,  at  some  point  It  Is  not  practical  to  continue 

reducing  the  voltage  as  the  signal  to  noise  ratio  will  become  Intolerable. 
Furthermore,  the  recently-discovered  funneling  phenomenon will  cause  the 
effective  charge-collection  volume  to  scale  differently  thm  the  pkgrslcal 
dimensions  of  the  volume  storing  the  charge  that  represents  a  bit.  Cbrrmt 
scaling  would  have  to  take  funneling  Into  account  as  the  device  dlmlnel^s 
become  smeller.  At  the  present  time  It  Is  not  clear  which  of  uwny  possible 
scaling  scenarios  should  be  used  for  these  exploratory  calculations.  Msdice, 
for  part  of  the  malysis  of  soft-upset  rates  reported  here,  a  wbretMutse 
assumption  has  been  made  that  aE  varies  as  Me  have  also  exi^bm  the 
effect  of  variation  of  scaling  laws  by  repeating  the  calculations  for  aE 
varying  as  a*^,  probably  the  most  realistic  scaling  assaalftloil  with  our 
present  knowledge. 

With  the  assumption  that  aE  scales  as  a'^,  and  the  scaling  of 
linear  dimensions  at  a~>,  the  minimum  chord  length,  tnin*  fw*  deposit  of 
energy  greater  than  aE,  scales  as 

(15) 

Furthermore,  the  lower  limit  of  Integration  Lq  scales  as 

1*0  •  Lo/«^»  ^11) 


Assuming  that  the  same  chip  ar«a  will  be  devoted  to  memory.  If 
the  reference  size  Is  M,  the  new  size  M*  Is 

M*  a  «2m.  (17) 

However,  the  quantity  S  in  equation  3  Is  the  surface  per  cell  so  that 

S*  -  S/(i2  (13) 

for  the  calculation  of  events  per  bltnJay. 

The  basic  equation  for  calculation  of  event-rate  after  scaling 

becomes. 


Ne  -  S/4a2  J  d{L)C(aE/«3pL)dL.  (19) 

•-0 

The  error  rate,  Ng,  Is  obtained  by  multiplication  by  c,  the 
error  conversion  factor. 

Ne  a  Ne  c.  (20) 

This  error  conversion  factor,  c,  which  depends  on  the  memory  configuration, 
is  (number  of  cells/memory  unit)  x  (vulnerability/cell). 

B.  EXAMPLES 

Application  of  the  preceding  theory  requires  availability  of  the 
cell  dimensions  and  critical  energy  for  specific  devices.  Following  Burke, 
we  use  the  parameters  developed  by  Pickel  and  Blandfordl3,18  for  a  number 
of  device  types.  The  parameters  used  in  the  calculations  are  listed  In 
Table  I. 

The  critical  dimensions  for  the  NMOS  dynamic  RAM  were  Inferred  by 
Pickel  and  Blandford  from  the  manufacturer's  data  on  the  device*  The  error 
conversion  factor,  c.  In  equation  20  Is  set  equal  to  1/2,  assuming  -  1/2  the 
cells  are  empty  at  one  time,  as  did  Pickel  and  Blandford. 

The  parameters  for  a  CMOS-Bulk  RAM  example  were  also  obtained 
from  reference  13.  Here  the  value  of  aE  -  22.5  MeV  Is  obtained  from 
heavy-ion  upset  measurements.  In  this  case,  since  there  are  six  devices  per 
memory  unit  and  1/2  are  vulnerable,  c  Is  set  equal  to  3. 

Similarly,  the  dimensions  for  the  4K,  CMOS-SOS,  five-transistor 
memory  cell  RAM’s  were  obtained  from  the  work  of  Pickel  and  Blandford. 
Following  the  approximation  used  by  Burke,  c  Is  set  equal  to  5  In  place  of 
the  averaging  procedure  used  by  Pickel  and  Blandford. 

The  results  of  calculation  of  soft-upset  rate  for  the  three  device 
types  In  the  90  percent  worst  case  spectrum  environment  are  presented  In 
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Figure  13.  As  In  all  of  our  calculations,  the  scaling  goes  from  a  -  0.5  to 
a  -  100,  well  beyond  the  VHSIC  region.  The  relative  vulnerability  of  the 
various  devices  Is  Indicated.  Furthermore  we  see  that  as  larger  scale 
Integration  proceeds,  the  predicted  soft-upset  rate  peaks  at  the  256K  to  1M 
memory  size  and  then  falls  off. 

The  dependence  of  predicted  soft-upset  rate  on  the  specific 
environment  Is  shown  In  Figures  14  to  16.  These  results  suggest  that 

detailed  analysis  of  the  environment  associated  with  a  particular  satellite 
orbit  Improves  the  predictive  accuracy  of  soft-upset  calculations. 

The  Heinrich  spectrum  used  by  Burke  In  his  calculations  Is  very  close 
to  the  solar  maximum  spectrum  used  In  the  NRL  calculations.  Although 
calculations  utilizing  the  solar  maximum  spectrum  yield  similar  results  In 
the  range  of  scaling  explored  by  Burke,  the  NRL  calculations  with  the  solar 
minimum  spectrum  agree  more  closely  with  Burke's  results.  The  fact  that 
there  Is  reasonable  agreement  between  the  Burke-  and  the  NRL-calculatlons 
Indicates  that  the  Burke  approximation  for  the  chord-length  distributions 

does  not  Introduce  any  large  errors. 

The  Importance  of  the  critical  energy  parameter,  aE,  Is  shown  In 
Figure  17.  Here,  the  critical  energy  has  been  arbitrarily  quadrupled  for  a 
calculation  of  soft  upset  rate  for  the  N-MOS  dynamic  RAM.  This  change  In  aE 
reduces  the  vulnerability  by  approximately  an  order  of  magnitude. 

The  calculations  presented  up  to  this  point  with  aE  proportional  to 

can  be  considered  a  worst-case  of  scaling  to  smaller  device  features. 

The  effect  of  scaling  the  critical  energy  as  a-2  Instead  of  o“3  is  shown 

In  Figures  18  to  20.  It  Is  quite  likely  that  the  «“2  scaling  Is  more 

applicable  to  the  way  device  scaling  will  go.  We  note  that  the  a-2 

scaling  predicts  soft-upset  rates  two  to  three  orders  of  magnitude  lower 
than  that  predicted  for  scaling  for  scale  factors  larger  than  10. 
This  strong  dependence  of  predicted  soft-upset  rate  on  scaling  scenario 
Indicates  that  a  more  detailed  Investigation  of  scaling  Is  required  for 

accurate  soft-upset  rate  predictions. 

The  predicted  soft-error  rate  due  to  direct  Ionization  by  protons  In 
the  proton  radiation  belt  at  600  nautical  miles  Is  sh^  In  Figures  21  and 
22  for  the  critical  energy  scaling  as  and  as  •*‘2.  These  results  are 
preliminary  as  work  Is  In  progress  on  the  soft-upset  rate  Induced  by  nuclear 
reactions  In  silicon. 

V.  SUMMARY  ANO  DISCUSSION 

Calculations  of  predicted  rates  of  soft-upset  failure  of  devices  In 
the  cosmic-ray  environment  are  presented  which  parallel  calculations 
performed  by  Burke.  The  present  calculations  utilize  improved  cosmic-ray 
environments  generated  at  NRL,  and  exact  calculation  of  the  Integral 
chord-length  distributions.  Furthermore,  the  scaling  Is  extended  to  smaller 
device  sizes. 

The  NRL  predictions  of  soft-upset  rate  yield  similar  results  to  those 


obtained  by  Burke.  This  agreement  indicates  that  the  Burke  approximation  to 
the  integral  chord-length  distribution  does  not  Introduce  large  errors. 

Calculations  utilizing  the  90  percent  worst,  case  spectrum  show  a  large 
increase  in  predicted  soft-upset  rate  in  the  scaling  range  of  approximately 
four  to  eighty  over  the  results  of  Burke.  The  90  percent  worst  case 
spectrum  is  to  be  preferred  for  a  more  realistic  estimate  of  the  soft-upset 
risk  for  a  satellite  in  the  cosmic-ray  environment  than  not  including  the 
contribution  due  to  low  level  solar  and  Interplanetary  activity.  Comparison 
of  soft-upset  rates  in  the  three  cosmic-ray  environments  utilized  in  these 
calculations  indicates  the  importance  of  accurate  evaluation  of  the 
environment  for  reliable  prediction  of  soft-upset  failure  rates. 

The  variation  of  predicted  soft-upset  rate  with  the  critical  energy 
shows  the  importance  of  correct  determination  of  the  device  parameters  for 
soft-upset  predictions.  A  calculation  with  two  different  scaling  scenarios 
shows  the  dependence  of  soft  upset  rate  on  the  details  of  scaling.  Clearly, 
further  investigation  of  scaling  is  required.  Perhaps,  it  would  be  more 
correct  to  use  different  scaling  scenarios  for  each  of  several  regions  of 
scaling. 

Preliminary  results  of  predicted  soft-upset  rates  in  the  proton 
radiation  belts  at  600  nautical  miles  are  presented.  These  calculations 
only  include  effects  due  to  direct  ionization  by  the  protons.  Errors  due  to 
direct  ionization  by  protons  in  the  radiation  belts  can  become  the  limiting 
factor  on  missions  as  devices  are  scaled  down.  For  the  examples  in  the 
present  calculations  there  is  a  rapid  increase  in  soft-upset  rate  at  a 
scaling  factor  of  approximately  four.  For  more  sensitive  devices,  failure 
rate  can  become  catastrophic  with  current  technologies. 

An  important  conclusion  that  can  be  inferred  from  these  exploratory 
calculations  with  several  scaling  scenarios  is  that  cosmic-ray  induced 
soft-upset  rates  do  not  increase  Indefinitely  as  feature  size  is  scaled  down. 
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Figure  1.  LET  spectrum  due  to  galactic  cosmic  rays  at  times  of  minimum 
solar  activity.  The  spikes  In  this  and  the  following  figures 
arise  from  singularities  In  the  conversion  from  energy  to  LET 
spectra.  It  Is  the  area  under  these  spikes  and  not  their  height 
that  carry  Information  on  the  differential  particle  flux. 
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Figure  2.  LET  spectrum  due  to  galactic  cosmic  rays  at  times  of  maximum 
solar  activity. 
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Figure  3.  90  percent  worst  case  spectrum.  LET  spectrum  due  to  galactic 

cosmic  rays  at  times  of  minimum  solar  activity  plus  90  percent 
worst  case  contribution  of  low  energy  particles  ou^o  solar  and 
interplanetary  activity. 
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Figure  4.  Cosmlc-rsy  LET  spectrtMi  behind  0  g/cm^  shielding  due  to 
Heinrich. 


Figure  11.  Integral  chord-length  distribution  for  3.5  x  14  x  21  microtaeter 
parallelepiped.  Comparison  of  exact  evaluation  and 
Approximation  2. 
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Figure  12.  Integral  chord-length  distribution  for  3  x  TO  x  10  aflcrOMeter 
parallelepiped.  Comparison  of  exact  evaluation  and 
Approximation  2. 
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Figure  14.  Prediction  of  soft-upset  rate  for  N-MOS  dynamic  RAM  in  three 
different  cosmic  ray  environments  for  aE  varying  as  o~3. 
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Figure  15.  Prediction  of  soft-upset  rate  for  CMOS  bulk  static  RAM  In  three 
different  cosmic  ray  environments  for  aE  varying  as  a"3. 
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Figure  16.  Prediction  of  soft-upset  rate  for  CMOS-SOS  static  RAM  In  three 
different  cosmic  ray  environments  for  aE  varying  as  a“3. 
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Figure  17.  Effect  of  change  In  critical  energy,  aE,  on  predicted  soft-upset 
vulnerability  of  N-MOS  d-RAM  In  90  percent  worst  case  spectrum 
environment.  Scaling  assumed  aE  varies  as  a~3. 
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Figure  18.  Effect  of  scaling  scenario,  aE  varying  as  a”3  vs  tE  varying  as 
o“2,  on  predicted  soft-upset  rate  for  N-MOS  dynamic  RAM. 
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Figure  19.  Effect  of  scaling  scenario,  aE  varying  as  oi-3  vs  aE  varying  as 
on  predicted  soft-upset  rate  for  CMOS  bulk  static  RAM, 
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Figure  21.  Prediction  of  soft  upset  rete  due  to  direct  lonlzitlon  by 
radiation  belt  protons  at  600  nautical  nlles  for  critical  energy 
scaling  as 
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Figure  22.  Prediction  of  soft  upset  rate  due  to  direct  Ionization  by 
radiation  belt  protons  at  600  nautical  miles  for  critical  energy 
scaling  as 
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TABLE  I.  Parameters  Used  In  Soft-Upset  Calculations 


Device  Type 

Device  Dimensions  for 

4k  Reference  Device 
(Scale  Factor  ■  1) 

(urn) 

Critical  Energy 
(MeV) 

Error  Conversion 
Factor 

N-MOS  d-RAM 

3.5  X  14  X  21 

5.6 

0.5 

N-MOS  d-RAM 

3.5  X  14  X  21 

22.5 

0.5 

CMOS-Bulk 
Static  RAM 

3  X  10  X  10 

22.5 

3.0 

CMOS-SOS 
Static  RAM 

0.5  X  5  X  15 

24.75 

5.0 
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